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ABSTRACT 

The velocity function derived from large scale surveys can be compared with the predictions 
of ACDM cosmology, by matching the measured rotation velocities Vrot of galaxies to the 
maximum circular velocity of dark matter (DM) halos Vmax- For Vrot<50kms“^, a major 
discrepancy arises between the observed and ACDM velocity functions. However, the manner 
in which different observational measures of Vrot are associated with Vmax is not straight 
forward in dwarf galaxies. We instead relate galaxies to DM halos using the empirical baryon- 
mass to halo-mass relation, and show that different observational measures of Vrot result in 
very different velocity functions. We show how the W50 velocity function, i.e. using the HI 
profile line width at 50% of peak HI flux to measure Vrot, can be reconciled with a ACDM 
cosmology. Our semi-empirical methodology allows us to determine the region of rotation 
curves that are probed by HI measurements (Rri), and shows that the Vrot of dwarfs are 
generally measured at a fraction of Rmax, explaining their tendency to have rising rotation 
curves. We provide fitting formulae for relating Rri and Reff (the effective radius) to the virial 
radius of DM halos. To continue to use velocity functions as a probe of ACDM cosmology, 
it is necessary to be precise about how the different measures of rotation velocity are probing 
the mass of the DM halos, dropping the assumption that any measure of rotational velocity 
can be equally used as a proxy for Vmax- 
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1 INTRODUCTION 

The rotation velocity Vrot of galaxies provides a measure of the 
mass within the radius at which the velocity measurement was 
made. The velocity function of galaxies, the number density of 
galaxies as a function of their velocity, is thus a close kin to the 
mass function, the number density of galaxies as a function of mass. 

The observed velocity function of galaxies has been derived 
by large scale HI surveys, using the HI line width as a measure 
of rotational velocity. The velocity (and mass) functions for dark 
matter halos in a ACDM cosmology can be readily derived for DM 
only simulations, with circular velocities defined as y/ GM(r)/r 
where G is the gravitational constant, and M(r) the mass inclosed 
within each radius r. Thus, one can compare the velocity function 
of observed galaxies to that of dark matter halos in simulations, by 
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matching the observed rotational velocities to the circular veloc¬ 
ity of dark matter halos. Generally, the maximum circular velocity, 
Vmax, is chosen as representative of DM halos, although the veloc¬ 
ity at any well defined radius could be chosen. 

The velocity function derived using the half width of the HI 
line profile, measured at 50% of the profile peak W50/2, has been 
shown to differ markedly from the Vmax function expect ed from 
ACDM cosmology for galaxies with W50/2 ^ 50kms~ ^ |Zavala 
et al.|2009[|Zwaan et al.|2010[|Trujillo-Gomez et al.|201 l[|Papaster- 
gis et al.|2011[|Klypinetal.|2014F 

However, for low mass galaxies, the relation between W50/2, 
and Vmax becomes less clear than for more massive disc galaxies 
for several reasons: low mass galaxies have a higher ratio of veloc¬ 
ity dispersion to rotation velocity, are thicker, with less regular HI 
discs, rising rotation curves which have not reached a maximum, 
and decreasing baryon fractions meaning that baryons may be con¬ 
fined to increasingly central regions of the DM halo. 

To make a more accurate comparison between observed and 
predicted velocities distributions |Papastergis et al.| ( |2015| ) relaxed 
the assumption that W50/2 and Vmax can be directly compared, 
and instead matched W50/2 velocities to the velocity of DM halos 
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Figure 1. The abundance matching relation adopted in this study, Mb- 
Mhaio (red line) and M*-Mhaio (black line). The region that is extrapolated 
from the Papastergis 2012 Mb-Mhaio relation is dashed. 
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Figure 2. The baryon mass versus W50/2 for ALFALFA galaxies (black 
dots) and the Local Volume catalogue of Karachentsev 2013 (red dots). 
These data sets share many galaxies. The data have similar distributions 
of low mass galaxies where the Local Volume statistics are better. 


at the largest radius to which the observed HI gas extends Riast • The 
mismatch between theory and observation was shown to persist, 
although this procedure has its own issues, such as the assumption 
of an NFW density profile for the DM halos, the assumption that 
concentration of the galaxies’ halos match those of pure N-body 
simulations, and the assumption that W50/2 refiects the rotational 
velocity of dwarf galaxies at Riast • 

In this study, we put aside the issue of how to best match ob¬ 
served and theoretical velocities, and instead match observed galax¬ 
ies to DM halos based on mass. Galaxies with regular, extended HI 
discs and well studied rotation curves, have been shown to follow a 
very tight relation between baryonic mass and circular velocity, the 
Baryonic Tully Fisher relation (BTFR). The small scatter in this 
relation provides compelling evidence that the baryonic mass of 
galaxies is tightly correlated with their total mass (e.g. |McGaugh| 
|2005| ). We propose that matching baryonic mass to halo mass pro¬ 
vides a more reliable mapping between observed galaxies and DM 
halos, than does the mapping between velocities. Several studies 
have explored the relation between galaxy masses, and the mass of 
DM halos (e.g. [Shankar et al.|200^ |Moster et al.|2010[ |Guo et al^ 
|2010[ [Papastergis et al.|2012| l 

Part of our motivation comes from the fact that the tight re¬ 
lation between baryon mass and rotation velocity when measured 
from the flat part of the rotation curve is not well reproduced, 
in terms of scatter, when one looks at the HI line-widths of low 
mass galaxies. For example, [Trachtemach et aL| ( |2009| ) explored 11 
dwarf galaxies with baryonic masses 0.14<Mb<12.92xlO^M©, 
and find that 6 have well defined rotation curves and do fit on the 
BTFR, while 5 others have W50/2 values that fall well below the 
BTFR: these 5 galaxies also have poorly defined rotation curves. 
When matching observations to DM halos by velocities, these lat¬ 
ter 5 galaxies would be placed in low mass halos. The study of 101 
dwarf galaxies in [Geha et ^ ( |2006| ), selected as low luminosity 
SDSS galaxies, also shows a large scatter in the values of W20 in a 
sample of galaxies with small range of baryon masses. 

The other point worth considering is that matching W50/2 to 
Vmax leads to many dwarf galaxies being matched to low mass ha¬ 
los, implying baryon fractions that are higher than what is found 
in L* galaxies. For example galaxies with W50/2 ~10kms“^, 
when assuming this reflects Vmax, are matched to halos with 


Mhaio~3x lO^M© (Figure 10 in Klypin et al.]2014) , even though 
they have Mb^lO^M©. 

Our contention is that the low value of W50/2 in many dwarfs 
is probing relatively inner regions of rising rotation curves, and 
therefore the Mb-Mhaio relation is a more robust manner for con¬ 
necting galaxies to DM halos. We will see in particular that dif¬ 
ferent measures of Vrot result in very different velocity functions, 
meaning that they cannot all be matched to Vmax equally. 

Matching galaxies to halos using mass provides an interest¬ 
ing, and we believe, compelling way to compare observed and the¬ 
oretical velocity functions. The study is organised as follows: in 
section we outline our procedure, then in section we show the 
adopted empirical relation between Mb and Mhaio and then show 
how Mb is related to Vrot in various observational data sets using 
various measures of Vrot- We then combine these relations, and 
show how differences in the Mb -Vrot relation affect the derived 
velocity function. We conclude that ACDM is compatible with the 
velocity function observations, under the condition that it can ex¬ 
plain the different BTFRs. In an appendix, we show that our results 
are robust to different data sets. 


2 METHODS 

Our procedure for deriving the velocity function is as follows: 

• use an empirical abundance matching between Mhaio and Mb. 

• use observed data sets to relate Mb to Vrot where the relation 
is dependent on the specific manner in which rotation velocity is 
measured (using HI line widths or rotation curves). 

• for each Mhaio we thus have empirically matched Mb and Vrot • 

• we then use the mass function from N-body simulations which 
gives us the predicted number of halos of each mass, and hence 
the predicted number of galaxies with particular values of Mb and 
Vrot, the latter being the velocity function. 

In our procedure, galaxies have been matched to halos according 
to their mass. We have not attempted to separate satellites and cen¬ 
tral galaxies which will be explored in a subsequent study; for the 
theoretical mass function we simply use a power law fit to the sim¬ 
ulations [Klypi££F^ ([20^ [20^. With galaxies matched to ha- 
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Figure 3. The baryonic Tully-Fisher relation (BTFR) for the Karachentsev 
2013 galaxies, with a fit to galaxies which have inclinations inc>60°. We 
use V 50 , the W50 line half widths adjusted for inclination and velocity 
dispersion. Also shown are the BTFRs of McGaugh 2012, derived using 
Vflat and from Gurovich 2010, derived using W20. 


los using the adopted Mhaio and Mb relation, the resultant velocity 
function differs depending on the adopted relation between Mb to 
Vrot i.e the BTFR, which itself differs according to the manner in 
which Vrot is measured. 

For our purposes, we are interested in adopting a single rela¬ 
tion between each of the galaxy properties Mhaio, M>, and Mgas 
and then exploring the effect of convolving these relations with 
different BTFRs. In the appendix we start with a different set of 
relations between Mhaio, M* and Mgas based on the recent study 
of [Bradford et al.| ( |2015| ), and then repeat our procedure. The re¬ 
sults are very similar to those presented in the paper, indicating that 
our results are robust. Within the paper, we also make several other 
consistency checks between different data sets. We neglect scatter 
in our adopted scaling relations as it will not affect our final deter¬ 
mination of the low-velocity end of the velocity function. 


3 ANALYSIS 
3.1 Mb-Mhaio 

We take a fit to the Mb-Mhaio relation from jPapastergis et al.j 
( |2012j ), using the equation: 


Mb = 0.0635Mhaio 


Mhaio (Mhaio 

1011.65 j + ) 


( 1 ) 

This has the same functional form as the M*-Mhaio relation as 
[Moster et al.j { |2010j ). The relation is shown as a red line in Fig¬ 
ure 1, with the region that is extrapolated from the jPapastergis et alT] 
\2012\ relation shown as a dashed line. We then compute the gas 
and stellar composition of the baryons using equation 1 of Papaster- 
gis 2012, 


logio(MHi/M.) = -0.431ogio(M.) + 3.75 (2) 


The adopted relation between gas and stellar mass of galaxies is 
also consistent with the empirical relation derived by [Peeples &| 
[Shankar [ ( [2011[ ). The residual M*-Mhaio relation, shown as a black 


Figure 4. HIPASS data showing the difference between HI line width mea¬ 
sured at 20% and 50% of the peak flux. We plot W50/2 versus W50AV20. 
For low mass galaxies, there is a big difference in the two indicators of ro¬ 
tation velocity, with a large tail toward low values of W50AV20. The fitted 
red line will be used to show the change in the velocity function when us¬ 
ing W20, as compared to W50, as an indicator of rotation velocity. The red 
dashed line shows the result of using W20 = W50+25 km s“^. 


line in Figure 1, is close to the [Guo et al. p010[ ) relation. It has 
been shown that the extrapolation of [Guo et al.[ ( [2010[ ) to lower 
mass galaxies gives a good approximation of abundance matching 
results for Local Group galaxies ( [Brook et al.[20I4[ l. 

Galaxies are matched down to halo masses of around 7x10®. 
Below this limit, several studies have shown that galaxy formation 
is suppressed by the early re-ionisation of the intergalactic medium 
( [Bullock et al.|200^[Somerville[2002[[Benson et al.|2002[ l. 


3.2 Mb-Vrot: measures of the BTFR 


We first explore the relation between W50/2 and baryonic mass for 
ALFALFA galaxies ( [Haynes et al.[201 ij), and for th e Local Volume 
(LV) galaxies as compiled by Karachentsev et aljj 2013[). W e cal¬ 
culate M* for the ALFALFA galaxies following |Bdl et al.[ ( [2003[ ) 
and calculate HI mass Mhi =2.36e5d^ xflux, where flux is the HI 
fiux and d is the distance in Mpc, then Mb=M*-H4/3MHi. We plot, 
in Figure 2, Mb versus W50/2, where the ALFALFA galaxies are 
shown as black circles and the LV galaxies as red circles. The dis¬ 
tributions are similar for the low velocity region we are interested 
in, and in fact there is significant overlap in the data sets. 

We define 


Vso = v/(W50/2/sm(inc))2 _ 32 


(3) 


taking account of inclination and velocity dispersion in deriving 
a measure of rotation velocity. In Figure 3 we then take the LV 
galaxies and plot the V 50 BTFR (red dots) and use galaxies that 
have inclinations inc>60° (red circles) to make a fit (black line): 


logioMb = 5.19 + 2.221ogioV5o (4) 

Also shown in Figure 3 are two BTFRs taken from the litera¬ 
ture; the BTFR derived using the fiat part of the rotation curve Vflat 
of galaxies with regular, extended HI discs and well resolved rota¬ 
tion curves, log^oMb = 2.01 + 3.821ogiQVflat ( [McGaugh|2012[ l; 
and the BTFR derived using inclination and dispersion adjusted 
W20 measurements, log^o^b = 3.7 + 31ogiQV2o ([Gurovich et al.[ 

[Mot . 
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Figure 5. The yellow line shows the Vmax velocity function of DM only 
simulations. We then show the velocity function predicted by ACDM 
cosmology when making different assumptions about the BTFR, having 
matched observed baryon masses to DM halo masses. Using the W50/2 
BTFR, equation!^ is shown as the red line. The W20/2 BTFR, derived us¬ 
ing equation|^is shown as a red dashed line. Results of assuming the BTFRs 
of Mcgaugh (2012) and Gurovich et al. (2010) are shown as blue solid and 
blue dashed lines. The observational results from Klypin et al. (2014) are 
shown as black circles with error bars 


3.3 W50 compared to W20 

Here, we highlight the difference between W50 and W20, in low 
mass galaxies, i.e. the difference between measuring the HI line- 
widths at 20% and 50% of the peak value. We use data taken from 
the HOPCAT catalogue ( [Doyle et al.|2005| ) of the HIPASS survey, 
but also show the relation from the recent study of [Bradford et al.j 
( [2015[ ), using the ALFALFA catalogue. 

Figure 1^ shows the ratio W50AV20 as a function of W50/2. 
For low mass galaxies, there is a significant difference in the two 
indicators of rotation velocity, with a large tail toward low values 
of W50AV20. We will use the fitted red line, 

W50/W20 = 0.1 + 0.41ogio(W50/2) (5) 

to adjust our velocity function from one based on W50, to one 
based on W20, shown as dashed red lines in figures 4, 6 & 7, high¬ 
lighting the effect of choosing W50 over W20 as an indicator of 
rotation velocity. This adjustment is almost identical to simply us¬ 
ing W20=W50-i-25kms“^(dashed red line in Figure]^, which is 
the typical difference between these measures of Vrot according to 
both [Bradford et al.[ ( [2015[ ) and [Koribalski et al.] ( [2004[ ). 


3.4 The Velocity Function 

In Figure]^ we explore the effect of different BTFRs on the veloc¬ 
ity function, when the BTFRs are used with our adopted Mb-Mhaio 
relation. In each case Vrot values are adjusted to line of sight ve¬ 
locities Vios by averaging over all inclinations to al low direct com- 
parison with the observed W50/2 velocity functior ~|Klypin et al. 
[2014[[Papastergis et al.[2015'] l. 


Figure 6. We plot Vrot as a function of the Vmax of the matched halo, un¬ 
der different assumptions for the BTFR after matching galaxies to halos on 
the basis of mass. We use 2 different empirical BTFRs from the literature, 
shown as blue lines. We also use the fit to W 50 from figure 2 (red line) and 
the adjustment from W50 to W20 as shown in Figure 4 (red dashed line). 


The yellow line in Figure shows the Vmax velocity func¬ 
tion of DM only simulations, while the the velocity function de¬ 
rived from the W50/2 BTFR, equation]^ is shown as the red line. 
The velocity function derived from the W20/2 BTFR using equa¬ 
tions |^&|^is shown as a red dashed line. Results of assuming the 
BTFRs of ^cGau^ ( [20 1 2[ ) , and [Gurovich et al.[ ( [2010[ ) are shown 
as blue solid and blue dashed lines respectively. The effect of dif¬ 
ferent BTFRs on the velocity function is significant. 

Consider the difference between the velocity function derived 
using the Vfiat BTFR, i.e. the blue solid line, as compared to the 
yellow line which shows the Vmax velocity function predicted from 
pure N-body simulations. The difference is large, yet hydrodynam- 
ical galaxy formation simulations within a ACDM cosmological 
context have been shown to simultaneously match both the Mb- 
Mhaio and the Vflat BTFR relation (Santos-Santos et al. submitted), 
which is sufficient to reconcile the lines. 

The red line, derived using the V 50 BTFR, shows the W50/2 
velocity function, which is very similar to the observed W50/2 ve¬ 
locity function (e.g. [Klypin et al.[20I4[ shown as black cycles with 
error bars in Figure pT 


3.5 Relating various measures of Vrot to Vmax 

Having matched galaxies to DM halos based on mass, we can ex¬ 
plore the implied links between various measures of Vrot for the 
observed galaxies, and the Vmax of matched DM halos. This is 
shown in Figurej^which shows the results for the different BTFRs. 
The Vrot taken from the two empirical BTFRs, as seen in earlier 
plots, are shown as a function of the Vmax of the matched DM ha¬ 
los (blue lines). Results using the V50 BTFR, i.e. equation are 
shown as the solid red line, and the adjustment from W50 to W20, 
equation]^ as the red dashed line. 


^ We have not included a dispersion correction as the observed veloc¬ 
ity functions use W50/2 directly. We derived the V 50 BTFR having used 
an inclination and dispersion corrected W50/2, because such relation fol¬ 
lows a single power law, with lower scatter, than the W50/2 BTFR. The 
re-inclusion of the velocity dispersion when deriving the velocity function 


is what causes the upturn in the W50/2 or W20/2 velocity functions at low 
velocities. The final result will not change if the velocity dispersion is main¬ 
tained throughout, so long as one selects an appropriate (but more compli¬ 
cated) relation between baryon mass and measured W50/2. 
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Figure 7. The red line is the Rmax of DM halos, plotted against the value of 
W50/2 for the matched galaxies. The blue circles show the radial extent of 
the HI discs (Rri) versus W50/2 of the FIGGs (Begum et al. 2008) sample. 
The blue line takes the fit of Rhi versus HI gas mass from Begum (2008), 
and converts from HI mass to W50/2 using equations[^&|^ 


Figure 8. The ratios Rni/Rvirand Reff/Rvir for observed galaxies and their 
matched halos. Rhi comes from Begum (2008) while Reff from Bradford et 
al. 2015. The black line and grey area at 0.01^0.005 x R 200 , is the relation 
found by Kravtsov (2013) for Rgff of more massive galaxies. 


3.6 The extent of low mass galaxies 


Finally, we highlight the difference between the extent of HI gas in 
galaxies, and the Rmax of their matched DM halos, where Rmax is 
the radius at which Vmax is reached. In FigureW50/2 is plotted 
as a function of Rhi, defined as the radius at which the HI gas 
density reached 19 atoms cm“^, for the FIGGs sample of galaxies 
( [Begum et al.|2008 ) shown as blue circles, a representative subset 
of the compiled Karachentsev et al.| ( |2013j ) data. 

The red line in Figure|7] shows the Rmax of DM halos, plotted 
against the value of W50/2 for the matched galaxies. The blue cir¬ 
cles show the radial extent of the HI discs (Rhi) versus W50/2 for 
the FIGGs sample. 

The blue line in Figure uses equation 1 of [Begum et al.j 
( [2008[) which relates HI mass to Rhi, along with our equa¬ 
tions]^ &E1 This is a good consistency check for the data sets. 

In Figure we plot the ratios Rni/RvirUnd Reff/Rvir for ob¬ 
served galaxies and their matched halos. Rhi comes from the 
FIGGS sample Begum (2008) while Reff from the sample of Brad¬ 
ford et al. 2015. The black line and grey area at 0.015 =f0.005x 
R 200 , is the relation found by [Kravts'o^ ( [2013] ) for Reff of more 
massive galaxies. A constant ratio of Reff/Rvir does not hold for 
lower mass galaxies, which are confined to more central regions 
within their host halos. The ratios are well fit power laws 


i^Sloi^Hl/Rvir — 7.65 -|- O.GlogioMhalo 


( 6 ) 


and 


I0R1 nHeff /H V 


4 DISCUSSION 

When matching W50/2 velocities of observed galaxies to Vmax of 
DM halos, one finds a significant excess in the number of DM halos 
with Vmax ^ 50km s“^. This mismatch occurs where halos are too 
massive to have had galaxy formation prevented by re-ionisation; 
from this perspective the observed velocity function indicates that 
there is a “too big to fail” problem ( [Zavala et al.[2009[[Klypin et~n?] 
[2014[[Papastergis et al.[2015| l. 


In this study we instead match observed galaxies to DM ha¬ 
los based on baryonic mass, and explored the impact of the vari¬ 
ous forms of the BTFR on the velocity function that is expected in 
ACDM cosmology. The BTFR derived from W50/2 values taken 
from large surveys of galaxies, results in a very different prediction 
for the velocity function than the one inferred from the BTFR de¬ 
rived using the Vffat values of dwarf galaxies that have extended, 
regular HI rotation fields. The question of whether the observed 
velocity function is consistent with ACDM cosmology may be re¬ 
duced to the question of whether the various forms of the BTFR 
can be reproduced within a ACDM cosmology. 

Yet it is has already been shown that hydrodynamical cos¬ 
mological galaxy formation simulations are able to simultaneously 
match the Vffat BTFR and Mb-Mhaio relation (Santos-Santos et al. 
2015 submitted, see also [Brook et al.[[2012| [Aumer et al.[[2013[ ), 
which is sufficient to reconcile the observed Vffat velocity function 
with ACDM predictions. Therefore, when considering galaxy for¬ 
mation processes, the predictions of ACDM are not embodied by 
the yellow line that is derived from purely DM only simulations: 
simulated galaxies within a ACDM context are compatible with 
the observed Vffat velocity function. 

Whether the W50/2 velocity function is compatible with 
ACDM cosmology becomes a matter of understanding how mea¬ 
sures of W50/2 in a large sample of galaxies relate to measures of 
Vffat in galaxies with extended HI discs and regular rotation curves. 

When matching galaxies to DM halos by mass, the problem 
lies with understanding observational measurements of Vrot rather 
than a shortfall in the number of galaxies that are “too big to fail”. 
This may require that the halos have DM cores, rather than cusps. 
Yet we already have significant evidence that many dwarf disc 
galaxies having slowly rising rotation curves, possibly indicative 
of cored DM profiles. 

Various proposals for explaining the existence of cores have 
been proposed, both within the ACDM paradigm (e.g. [Ni 


et al.[[ 19961 [Read et ^[2006 [ [Governato et al.[[20T0l [Pontzen & 


Governato [2012[l, or by taking an alternative form of DM (e.g. 
Zavala et al.[2009[[Vogelsberger et al. [2012[ [Schneider et al.[2014[ |. 

In an accompanying paper (Brook & Di Cintio 2015 submitted), 
we will compare how galaxy cusps and cores are refiected in the 
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velocity function, and then proceed to highlight how different core 
formation mechanisms may leave signatures in the detailed shape 
of the velocity function. 

According to our analysis, Vrot is not measured at Rmax in 
dwarf galaxies. We provide relations between Rhi and Res of ob¬ 
served galaxies, and Rvir of their matched DM halos, and showed 
that dwarfs are confined to more central regions of their matched 
DM halos than more massive galaxies. This explains why rotation 
curves of dwarf galaxies are generally rising at the last measured 
point. As one does not reach Rmax the measured Vrot is dependent 
on the extent of the HI gas, as well as its detailed structure. Scat¬ 
ter begins to emerge in the various BTFRs for low mass galaxies, 
and the different measures of Vrot become increasingly different, 
leading to different velocity functions. 

We performed several consistency checks between different 
data sets from the literature. In the appendix, we have taken this 
further, redoing our analysis using the recent data from [Bradford] 
jeTaLl l jMS] ), who measured the HI line-widths of a large sample 
of dwarf galaxies selected from the SDSS catalog. The results, and 
in particular the large effect of different measures of Vrot on the 
velocity function, are unchanged from the analysis in the paper. 

Our study highlights that if the velocity function is to be used 
as a probe of cosmological paradigm, we need to be precise in relat¬ 
ing measures of Vrot with the Vmax of DM halos. Not all measures 
of Vrot are equal and not all BTFRs have the same slope, particu¬ 
larly in low max galaxies. In a forthcoming paper, we will explore 
various forms of the BTFR within hydrodynamical, cosmological 
simulations, looking at how the different measures of the BTFR can 
be reconciled within a ACDM context. 
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APPENDIX A: DIEEERENT DATA SETS 

We have emphasised in our paper the effect of the different mea¬ 
sures of the BTFR on the velocity function. Our analysis adopts 
a set of empirical relations between Mhaio, M*, Mgas and Vrot, 
using data from several different studies. Here, we repeat our anal¬ 
ysis but starting with a different set of empirical relations between 
Mhaio, M*, M^as and Vrot, using the results of the recent study 
of [Bradford et al.|f2015| ). who study HI line-widths of a sample of 
isolated dwarf galaxies selected from the Sloan Digital Sky Survey. 
We show that our main conclusions are not altered. 


Al Different baryonic abundance matching 

Our starting point is the empirical M>,-Mhaio relation of [Moste^ 
[et al.[ ( |20T0] ), shown as a black line in figure [Al[ wh ich is slightly 
fiatter at the low mass end than the [Guo et al.| ( |2010] ) relation used 
in the paper. We then use the double power law relation between 
M* and Mgas from [Bradford et al.] ( [2015| ), which results in the Mb- 
Mhaio relation shown as the red line in figure [Al| 
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Figure Al. The abundance matching relation adopted in the appendix, M^- Figure A2. The baryonic Tully-Fisher relation (BTFR) for the data set, with 
Mhaio (red line) and M*-Mhaio (black line). a fit to galaxies. Also shown are the BTFRs of McGaugh 2012, derived 

using Vflat- 


A2 Reconciling BTFRs from a new data set 

[Bradford et al.| ( |20l3] ) derive a BTFR using W20/2, 

logioW2o/2 = -0.668 + 0.2771ogioMb (Al) 


where there is no correction made for dispersion. This BTFR is 
shown as a black line in Figure |A2| and is is very similar to the 
BTFR from |McGaugh| ( [20T^ , shown as a blue line. We also show, 
as a red line, the V 50 BTFR from this paper, i.e. equation]^ The 
dashed red line adjusts this V 50 BTFR to a W20/2 BTFR using 
equation after also adjusting V 50 to W50/2 by including the 
8 km s“^ velocity dispersion. The resultant relation is very close to 
the V 20 relation of jBradford et al.| ( |2015| ). 

[Bradford et al.| ( [20T5] t find a typical difference bet ween W50 
and W20 is ~25km s“^(see also Koribalski et al.|20()^ . We there¬ 
fore use this value to also transform the [Bradford et al.[ ( [2015[ ) 
BTFR (black line) to a W50 BTFR, shown as the dashed black 
line. The adjusted BTFR is very similar to our adopted W50 BTFR, 
equation|^ This gives us significant confidence in our analysis, and 
in particular in equations equation [^&[^from our paper, which are 
the key results which are driving the significant differences in the 
velocity function for different measures of Vrot • 


A3 The effect of BTFRs on the Velocity Function 

Starting this time from the abundance matching by mass shown in 
Figure [ aT] we show in Figure [A3 [ the effect of the different BTFRs 
on the velocity function. The yellow line shows the Vmax velocity 
function of DM only simulations. Using the W50/2 BTFR, equa¬ 
tion]^ is shown as the red line. The W20/2 BTFR, derived using 
equation]^ is shown as a red dashed line. Results of assuming the 
BTFR of Bradford (2015), i.e. equation [ aT] are shown as the black 
solid line, with the adjustment of W50=W20-25 km s“^shown as 
the black dashed line. 



Figure A3. The yellow line shows the Vmax velocity function of DM 
only simulations. We then show the velocity function predicted by ACDM 
cosmology when making different assumptions about the BTFR, having 
matched observed baryon masses to DM halo masses. Using the W50/2 
BTFR, equation]^ is shown as the red line. The W20/2 BTFR, derived us¬ 
ing equation[^is shown as a red dashed line. Results of assuming the BTFR 
of Bradford (2015), i.e. equation [a7] are shown as the black solid line, with 
the adjustment of W50=W20-25 km s“^shown as the black dashed lines. 
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